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INTRODUCTION 

The  NAVAIRDEVCEN  (Naval  Air  Development  Center)  is  a  large  user 
of  vehicle  and  propulsion  design  and  performance  computer  codes,  A 
significant  number  of  these  programs  is  dependent  on  the  use  of  input 
tabular  data  sets.  Experience  has  shown  that  it  is  desirable  to  have 
a  visual  representation  of  these  data  sets  prior  to  use  in  these  programs 
to  avoid  employing  incorrect  data.  In  addition,  the  programs  often 
require  a  means  of  expressing  output  results  in  a  form  other  than 
digital  data.  This  present  report  describes  a  general  purpose  routine 
which  will  generate  data  plots  in  terms  of  one  or  two  independent 
variables.  This  code,  entitled  GPPR  (General  Purpose  Plotting  Routine) 
for  plotting  data,  was  developed  for  the  NAVAIRDEVCEN  CDC  6600/CYBER  175 
computer  facilities  and  is  used  in  conjunction  with  a  model  735  CALCOMP 
pen  plotter.  A  user’s  guide  for  this  code  is  shown  in  Appendix  A. 

A  sample  problem  is  illustrated  in  Appendix  B  and  a  FORTRAN  code  listing 
is  given  in  Appendix  C. 


DISCUSSION 

CODE  DESCRIPTION 

The  GPPR  code  was  developed  as  a  computer  tool  to  permit  users  to 
easily  and  quickly  plot  digital  data  for  both  use  in  reports  and  editing 
of  table  data  sets  used  as  inputs  in  other  computer  codes.  A  completely 
general  plotting  routine  was  thought  to  be  awkward  and  inconvenient 
to  the  user,  because  of  the  large  numbers  of  inputs  that  would  be 
required.  For  this  reason  decisions  were  made  which  restrict  the 
options  related  to  the  plot  size  and  axis  labeling.  The  code  package 
consists  of  subroutines  GPPR,  AXSCALE,  and  function  SPLNQ1.  A  user 
written  main  program  is  required  to  use  the  plot  package.  An  example 
of  a  main  program  is  shown  in  Appendix  C. 

The  basic  features  and  options  selected  for  the  GPPR  subroutine  are 
as  follows: 

Plot  Si2e 

The  standard  size  GPPR  plot  is  8.5  by  li  inches  which  is  compatible 
with  present  NAVAIRDEVCEN  report  page  size.  The  manner  in  which 
these  plots  are  produced  on  the  CALCOMP  Pen  Plotter  is  shown  in  Figure  1. 
The  entire  plot  size  (all  symbols,  scales,  and  other  characters)  may 
be  changed  from  the  standard  size  through  aninput  size  factor.  For 
example  plots  of  4. 25  by  5.5  inches  would  be  obtained  with  an  input  size 
factor  equal  to  .5. 
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FIGURE  1.  GPPR  OUTPUT  FROM  CALCOMP  PEN  PLOTTER 
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Titles,  Axis  Labeling  and  Increment  Size 

The  standard  size  axis  lengths  are  5  and  7  inches,  respectively,  with  tic 
marks  provided  at  one  inch  intervals.  The  scale  and  increment  size  for 
each  axis  is  determined  from  another  subroutine,  AXSCALE.  The  AXSCALE 
routine  determines  the  largest  scale  and  increment  size  that  will  permit 
all  data  points  to  lie  within  the  boundary,  .01  inch  less  than  the  start 
of  and  .05  inch  greater  than  the  end  of  each  axis.  The  AXSCALE  routine 
will  select  the  correct  scale  and  increment  size  in  multiples  of  1,  2, 
or  5.  This  particular  scale  selection  procedure  enables  the  user  to 
overlay  the  GPPR  plot  with  10  divisions  per  inch  graph  paper  and  accurately 
read  values  from  the  GPPR  plot.  For  example,  if  y  axis  values  ranged 
from  0  to  20(1  the  axis  label  values  would  be  50  units  per  inch  with  a 
full  scale  value  of  350.  A  FORTRAN  listing  of  the  AXSCALE  code  is  found 
in  Appendix  €• 

Each  axis  title  label  is  input  by  the  user.  A  title  label  option 
for  a  second  independent  variable  is  also  provided.  Each  line  of  the 
main  title  heading  is  automatically  centered  over  the  plot.  The  entire 
main  title  heading  may  contain  a  maximum  of  four  lines  .with  no  more  than 
35  characters  per  line. 

Data  Presentation 


The  GPPR  routine  has  three  options  concerning  presentation  of  data 
on  each  plot. 

1.  Option  1 

The  first  option  is  to  plot  symbols  for  each  input  data  point.  If 
two  or  more  second  independent  variables  are  input  then  different  type 
symbols  will  be  used  for  each  group  of  points  representing  a  constant 
value  of  the  second  independent  variable. 

2.  Option  2 

The  second  option  is  an  extension  of  the  first  option.  The  symbols 
are  plotted  as  in  option  1  and  an  in-house  developed  cubic  spline  inter¬ 
polation  code,  function  SPLNQ1,  based  on  methods  in  reference  (a),  is 
used  to  draw  a  smooth  curve  between  each  set  of  symbol  types.  The  cubic 
spline  interpolation  technique  is  unique  in  that  the  interpolated  curve 
passes  through  each  data  point  and  has  continuous  first  and  second 
derivatives.  A  FORTRAN  code  listing  for  function  SPLNQ1  is  given  in 
Appendix  C. 

3.  Option  3 

The  third  option  permits  the  user  to  plot  the  interpolated  spline 
curve  with  symbols  appearing  only  at  each  end  of  the  spline  curve.  The 
symbols  at  the  end  of  each  spline  curve  are  retained  only  to  identify 
eaeh  curve  in  terms  of  the  associated  second  independent  variable  value. 
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4.  Grid  Option 

The  grid  option  permits  a  1  inch  grid  to  be  drawn  at  the  tic  marks  on 
the  axes.  Instructions  for  using  these  options  are  given  in  Appendix  A. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  general  purpose  plotting  routine  discussed  in  this  report  is  a 
valuable  computer  based  tool.  The  GPPR  subroutine  is  flexible  and 
easily  incorporated  into  new  or  existing  computer  codes. 

The  present  structure  of  GPPR  restricts  usage  to  curves  with  open 
arcs  and  single  values  of  the  dependent  variable  at  each  independent 
variable.  It  is  recommended  that  program  development  proceed  which  will 
extend  GPPR  to  permit  the  plotting  of  data  which  can  be  the  form  of  a 
closed  arc. 


REFERENCES 

(a)  Pennington,  Ralph  H. ,  "Introductory  Computer  Methods  and  Numerical 
Analysis",  Macmillian  Company,  London,  1970 
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APPENDIX  A 
USER'S  GUIDE 


All  inputs  into  the  GPPR  routine  are  passed  as  formal  parameters  in 
the  call  statement  to  GPPR.  The  order  of  the  parameters  in  the  call 
statement  is  as  follows: 

CALL  GPPR  ( NPLOT ,  LABY,  Nl,  LABX,  N2,  LABTL,  NT,  X,  Y,  NPTS,  LABVAL, 

NCC,  VLABL,  NDECV ,  ITIP,  IGRID ,  FAC) 

These  parameters  are  defined  as  follows: 

NPLOT  is  an  initializing  parameter  which  is  set  to  0  before  the 
first  call  to  GGPR.  NPLOT  is  only  set  once. 

LABY  is  an  array  containing  the  Hollerith  data  for  the  Y  axis  label. 

"Nl  is  the  number  of  elements  in  the  LABY  array  comprising  the  Y  axis 
label.  (1  element  **  10  characters) 

LABX  is  an  array  containing  the  Hollerith  data  for  the  X  axis  label. 

N2  is  the  number  of  elements  in  the  LABX  array  comparisons  for  the  X 
axis  label  (1  element  =  10  characters).  , 

LABTL  is  an  array  containing  the  Hollerith  data  for  the  main  title 
labeTi  Groups  of  3  or  more  consecutive  blank  characters  between  non 
blank  characters  indicate  the  beginning  of  a  new  line  in  the  title 
label.  Space  limitations  permit  a  maximum  of  4  lines  with  approximately 
35  characters  per  title  line.  If  the  user  attempts  to  use  a  title  line 
longer  than  35  characters,  the  line  will  terminate  at  the  next  blank 
after  the  thirty-fifth  character  in  that  line.  Leading  and  trailing 
blank  characters  for  each  line  are  ignored  and  each  line  of  title  is 
centered  on  the  plot. 

NT  is  the  number  of  elements  in  the  LABTL  array  comprising  the  main 
title  label  (all  lines)  (1  element  *  10  characters).  If  NT  is  equal  to 
0,  no  main  title  is  written. 

X  is  an  array  containing  the  X  axis  values  of  data  for  the  entire 
plotT  For  example,  consider  a  plot  with  three  curves  to  be  drawn  as 
shown  in  Figure  A-l.  If  the  first  curve  has  6  data  points  and  the 
second  curve  has  4  data  points,  then  elements  X(l)  to  X(6)  are  the  X 
axis  coordinate  values  for  the  first  curve  and  elements  X(7)  to  X(1Q) 
are  the  X  axis  coordinate  values  of  the  second  curve.  The  third  curve 
points  would  follow  in  X(ll)  to  X(15).  The  order  of  the  6  points  in  the 
X  array  comprising  the  first  curve  is  unimportant.  Similarly,  the  order 
of  the  points  comprising  any  curve  is  unimportant. 
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Y  is  an  array  containing  the  Y  axis  values  of  data  for  the  entire 
plot.  There  must  be  a  one  to  one  correspondence  between  the  Y  array 
values  and  the  X  array  values.  For  example,  X(l)  and  Y(l)  are  the  X  and 
Y  axis  coordinate  values  of  the  first  input  data  point. 

NPTS  is  an  array  containing  the  number  of  data  points  for  each  curve 
on  the  plot.  In  the  preceeding  example,  NPTS  (1)  =  6,  NPTS  (2)  =  4,  and 
NPTS  (3)  =  5.  The  NPTS  element  following  the  last  element  used  to 
identify  the  last  curve  must  be  set  equal  to  zero.  In  this  example, 

NPTS  (4)  must  be  set  equal  to  0. 

LABVAL  is  an  array  containing  the  Hollerith  data  for  the  second 
independent  title  label.  This  title  label  is  used  identify  the  value 
held  constant  for  each  curve  on  the  plot.  (1  element  of  the  array  will 
contain  10  eharacters.j 

NCC  is  the  number  of  elements  in  the  LABVAL  array  comprising  the 
second  independent  variable  (if  NCC  «  0  no  second  independent  variable 
title  is  generated  on  the  plot)  title  label  (1  element  =  10  characters). 

VLABL  is  an  array  containing  the  values  of  the  second  independent 
variable  associated  with  each  curve.  The  first  element  of  VLABL  is 
associated  with  the  first  set  of  points  in  the  X  and  Y  arrays.  If  NCC 
is  set  equal  to  zero  then  VLABL  is  not  applicable. 

NDECV  is  the  number  of  significant  figures  to  the  right  of  the 
decimal  for  the  VLABL  array  to  be  used  on  the  plot. 

ITIP  is  a  switch  used  to  indicate  the  following  options. 

ITIP  OPTION 

1  plot  symbols  only 

2  plot  symbols  and  draw  a  spline  curve  fitted 

with  respect  to  X  axis 

3  same  as  ITIP  *  2  except  symbols  are  only 
plotted  at  end  points  of  spline  curve 

-2  plot  symbols  and  draw  a  spline  curve 

fitted  with  respect  to  Y  axis 

-3  same  as  ITIP  *  -2  except  symbols  are  only 

plotted  at  end  points  of  the  spline  curve 

IGRID  is  a  switch  used  to  indicate  grid  options.  If  IGRID  =  1»  a  one 
inch  grid  is  drawn  on  the  plot.  Otherwise,  no  grid  is  drawn. 
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FAC  is  a  value  indicating  the  size  factor  of  the  plots.  FAC  set  to 
1.0  is  the  standard  size  indicating  plots  8.5  by  11  inches  will  be  drawn 

The  last  step  in  any  program  using  the  GPPR  subroutine  is  to  end  the 
plotting  tape.  To  end  the  plotting  tape  the  user  must  enter  a  call  to 
GPPREND  using  the  same  formal  parameters  previously  defined  in  the  call 
to  the  GPPR  subroutine. 
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APPENDIX  B 
SAMPLE  PROBLEM 


In  this  appendix  an  example  problem  using  the  GPPR  subroutine  is 
illustrated.  In  this  example,  a  2  independent  variable  function  is 
plotted,  standard  size.  In  order  for  the  user  to  implement  the  plotting 
package,  data  must  be  input  into  the  GPPR  subroutine  from  another  program 
or  subroutine.  In  this  example,  program  EXPLOT  is  the  means  by  which 
the  plot  inputs  are  transferred  into  the  GPPR  subroutine.  In  Appendix  C 
the  FORTRAN  listings  of  program  EXPLOT  and  the  other  routines  forming 
the  plotting  package  (GPPR,  AXSCALE,  and  SPLNQ1)  are  shown.  On  cards  3 
and  4  of  EXPLOT  (see  page  C-2),  the  input  data  points  (X  is  the  first 
independent  variable,  Y  is  the  dependent  variable)  are  set  using  a  DATA 
card  statement.  The  NPTS  array  is  set  to  7,  5,  and  0  on  card  6  indicating 

two  curves  are  input.  The  first  curve  has  7  points  (X  and  Y  array _ 

elements  1-7)  and  the  second  curve  has  5  points  (X  and  Y  array  elements  8- 
12).  The  elements  of  the  second  independent  variable  array  are  set  on 
card  5,  page  C-2,  to  values  ,5  and  .821  respectively.  The  axis  labels, 
main  title  label  and  second  independent  variable  label  are  set  on  cards 
7  through  11,  page  C-2.  The  call  to  the  GPPR  subroutine  is  shown  on 
card  12.  Most  of  the  formal  parameters  have  been  defined  above.  The 
remaining  parameters  are  defined  with  numerical  values  in  the  call 
statement  to  GPPR  indicating  the  following: 

The  Y  axis  label  has  2  elements 

The  X  axis  label  has  2  elements 

The  main  title  label  has  6  elements 

The  second  independent  variable  label  has  3  elements 

The  values  of  the  Z  array  on  the  plot  will  have  3  significant 
figures  to  the  right  of  the  decimal 

The  plot  will  be  standard  size  with  symbols,  spline  curve,  and 
grid  drawn 

The  resulting  plot  generated  from  the  CALCOMP  Pen  Plotter  is  shown 
in  Figure  B-l.  The  computer  time  required  to  generate  the  CALCOMP  Pen 
Plotter  instructions  for  this  problem  was  about  .6  CP  (Central  Processor) 
seconds. 
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PROGRAM  EXPLOT 
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